Cooling of photoexcited carriers in graphene by internal and substrate phonons 
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We investigate the energy relaxation of hot carriers produced by photoexcitation of graphene 
through coupling to both intrinsic and remote (substrate) surface polar phonons using the Boltzmann 
equation approach. We find that the energy relaxation of hot photocarriers in graphene on commonly 
used polar substrates, under most conditions, is dominated by remote surface polar phonons. We also 
calculate key characteristics of the energy relaxation process, such as the transient cooling time and 
steady state carrier temperatures and photocarriers densities, which determine the thermoelectric 
and photovoltaic photoresponse, respectively. Substrate engineering can be a promising route to 
efficient optoelectronic devices driven by hot carrier dynamics. 



I. INTRODUCTION 



Upon fast excitation of graphene carriers with light or other means, the dynamics of the resulting non-equilibrium 
carrier distribution evolve on a fast time scale and has been extensively studied both experimentally [THE] and 
theoretically frjJflS]- The relaxation involves an initial fast evolution towards quasi-thermal distribution on a 
femtosecond timescale via electron-electron collisions [TH EE ED HE], followed by energy transfer to phonons on a 
longer picosecond timescale. The conversion of the excess energy of these photoexcited carriers into electrical current 
before they lose this energy to the phonon baths represents one of the key challenges to efficient optoelectronic device. 

In this paper, we study the energy relaxation pathways of the photoexcited carriers via different inelastic scattering 
channels. Energy relaxation processes in graphene due to intrinsic optical and acoustic phonons have already been 
studiedEil [T^I21j . High energy optical phonon emission by hot carriers is responsible for the subpicosecond fast 
cooling process[TJ 121 122], followed by cooling via the acoustic modes. The latter is a slow process, that creates an 
electron-phonon cooling bottleneck|19j. Here, we focus on an extrinsic mechanism for cooling of photoexcited carriers 
in graphene via the remote surface polar phonon modes (SPP) of the substrate and compare their efficiency under 
different conditions with those of the internal phonon modes. 

In polar substrates such as Si02, a non- vanishing fluctuating electric field is generated by the propagating surface 
phonon modes [23 . The interactions of these SPP modes with charged carriers in the conduction channel was first 
explored in the context of inversion layer of semiconductor-oxide interface [24 26 . They have also been studied in 
other material systems such as carbon nanotubes |27If2"5] . where close proximity between charged carriers and the 
underlying substrate renders the SPP-phonon scattering more prominent. Similarly, in graphene, SPP was found to 
limit electronic transport properties [30 -39 . Recently, the SPP coupling with graphene plasmons was also probed 
experimentally through infrared spectroscopy[40-42 . In this work, we found irrespectively of the mechanism i.e. 
thermoelectric or photovoltaic, that SPP limits the overall strength of the steady state photoresponse on common 
substrates, and our results suggest that elimination of the SPP cooling channel can lead to an order of magnitude 
enhancement in the photoresponse. 

In Sec.[Hj we present the general theory, where details of the models for the electron cooling power via the different 
phonon baths are presented in the Appendix. We present the results of our calculation of the cooling powers in 
Sec. Ill A and discuss their relative contribution in detail, as a function of doping and electronic/lattice temperatures. 
In Sec. HTjB, we apply the above models to the study of the cooling dynamics of hot carriers due to continuous or 
pulsed light excitations. We calculate key experimental observables such as the transient cooling time, and steady 
state quantities such as the non-equilibrium electronic temperatures, excess photocarriers density and the out-of-plane 
thermal conductivity for graphene on common substrates. 
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II. THEORY AND MODELS 



Transition probability for emission and absorption of phonons with a particular phonon bath a is described by the 
Fermi's golden rule, 

S Q (k,k') = ^ l M k,k<| 2 {^A'-k-qW -Ek~ faq) + {N u , + l)^- k+q 5(K - £ k + hw q )} (1) 

q 

where q is the phonon momentum, N u = [exp(huj q / 7c_b7l) — is the Bose-Einstein distribution and Mg k , are the 
transition matrix elements related to the coupling with phonon bath a, to be defined below. For brevity, summation 
^ q <5k'-k±q shall be implicit hereafter. The cooling power is computed numerically by accounting for the transfer 
of electronic energy to the lattice during each scattering event i.e. E^ — E^. Therefore the net cooling power is 
calculated viaOH ED S3 M , 



K S ' = ^EE^( k ' k ')(^-SkO/k(l-/kO 

k k' 

= ^7§? E l M k,k'| 2 <5(£ k - - nw q )(E k , E k )F(k, k') 

k,k' 

pQO />00 /-27T 

= / gsg " / kdk / fc'dfc' / d6»|Af^ k ,| 2 5(fc' - k - UqvZ 1 )^ 1 - k)F(k,k') (2) 
(2Tr) 2 hJ J J Q ' 

where 2l and Tg are the lattice and electron temperatures respectively, and the electron distribution function is 
described by /k, and we define a composite Fermi-Boson distribution function as, 

J-(M') = (JVo,, + l)/k'(l - h) - N u J k (l - / k (3) 

As indicated by experiments J] [2] [18| , the electronic system is thermalized by the electron-electron interactions which 
occur at much faster timescale than the electron-phonon processes we are calculating here. Hence, it is appropriate 
to simply assume that /k follows the Fermi Dirac distribution function, i.e. [1 + exp(/3(£'k — where f3 — l/fesTg 

and /i is the chemical potential, controlled by chemical doping or electrical gating in experiments. It is apparent that 
the composite electron-phonon distribution function F(k, k') becomes zero when Tg = Tl, hence zero cooling power. 



In this work, we are interested in the energy exchange of electrons with the different phonon baths i.e. intrinsic 
acoustic phonons (AP), optical phonons (OP) and surface phonon polaritons (SPP). Vibrations of the substrate ions 
with the opposite charge polarity produce an electric field which decays exponentially away from the surface. Carriers 
in the nearby graphene can feel this electric field and be scattered by the SPP phonon. The decay length of the electric 
field is determined by the momentum transfer in the electron-phonon scattering event. For typical carrier density in 
graphene, the relevant momentum transfer is of the order of nm _1 , such that a substantial coupling strength is expected 
for graphene placed at van der Waals distance of 3.4 A away from the substrate. The transition matrix elements, 
Mkk', for electron interaction with AP, OP and SPP phonons are well-known in the literatures [26, 32, 33"! I45ti47] , 
We therefore defer their discussions to the Appendix, with the parameter set summarized in Table 1. We shall focus 
on the key results in what follows. 



III. RESULTS AND DISCUSSIONS 



A. Competing cooling pathways 



We begin with a simple illustration of the possible cooling pathways for photoexcited carriers in graphene in 
Fig.[T] Each thermal bath can be characterized by their respective temperatures T Q , and are in general different 
from the ambient temperature To. The heat exchange between these thermal baths can be described by the thermal 
conductivity, k, defined as the ratio between the power exchange per unit temperature difference i.e. 5P/5T. In 
general, the different phonon baths can each establish a different temperature upon interactions with the electrons 
(see also discussion in Sec. Ill B I. In this work, we shall assume a common temperature for all these phonon baths, 
denoted simply as the lattice temperature Tl. 



For a typical Si02 substrate thickness of ft, = 50 — 300 nm, k = n/h varies in a range w 5 MW/Km 2 — 10 MW/Km 2 , 
where Si02 film thermal conductivity is k = 0.5 — 1.4W/mK 49J. The interface thermal conductance of graphene 
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TABLE I: Parameters for the optical (OP), acoustic (AP) and substrate phonons (SPP). For SPP, we consider SiC>2 and h-BN 
substrates. ei ow (thigh) is the low (high) frequency dielectric constant of the dielectric and the surface optical phonon (SO) 
energies are obtained from the bulk longitudinal optical phonon (LO) phonons as hujso = ^-LoC txt/^ ) 1//2 - In this work, 
we consider only the two SO modes with the strongest coupling strength, denoted as uj\ and UJ2- Fj is the electron coupling 
parameter with the SO modes. For internal phonon modes, the energies of OP (r and K) and their deformation potential D op 
used in this work are summarized, and the sound velocity vs and deformation potential for the AP mode, D ac . See also the 
respective Appendix for details. 



SPP 


Si0 2 |33] 


h-BNg8] 


OP 


r 




K 




3.9 


5.09 


D op (eVA" 1 ) 


11 




16 


t-high 


2.4 


4.1 


fkuo (meV) 


197 




157 


Tujj\ (meV) 


58.9 


101.7 










hoj2 (meV) 


156.4 


195.7 




AP 






Fl (meV) 


0.237 


0.258 


D ac (eV) 




7.1 




Fi (meV) 


1.612 


0.52 


v s (km/s) 
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on SiC>2 substrate has been measured using various experimental techniques [50- 54], with values ranging from 
25MW/Km 2 — 180MW/Km 2 . On the theory front, several approaches have been employed to estimate this 
interface thermal conductance ,55 -_59 : , which varies from « 1MW/Km 2 - 100MW/Km 2 . As illustrated in Fig.jlJ 
energy transfered from electrons to the internal phonon baths is conducted to the underlying substrate through a 
phonon- limited ktb- i^tb between carbon surface and SiC>2 substrate has been estimated from molecular dynamics 
and is « 60 MW/Km 2 [55], and can depend also on the surface roughness. Alternatively, energy can be transfered 
directly to the substrate via the SPP phonons, i.e. kspp, and can depend sensitively on doping. For an undoped 
graphene, k-spp is on the order of 1MW/Km 2 while klat is even smaller, as we will see later in the discussion. In 
this section, we discuss how these cooling pathways depend on the various experimental conditions. 

Detailed balance condition of in- and out-scattering processes requires that V a vanishes under equilibrium 
condition i.e. ST = Te — Tl = 0. In the theory, this is ensured by the composite Fermi-Boson distribution function 
F(Te,Tl). The energy exchange efficiency with these various phonon baths depends upon, among other factors, the 
doping and electronic/lattice temperatures. Using the models described in Sec.|n| we calculate V a (TE) due to the 
various phonon baths for intrinsic/doped graphene under cold/hot (defined at Tj, = 10, 300 K respectively) lattice 
temperature as shown in Fig. [2] 

First, we discuss results on intrinsic graphene, see Fig.[2]^a,c,e), which can be understood on the basis of scattering 
phase space arguments. For cold neutral graphene, Pauli blocking limits the electronic transitions involved to mainly 
interband processes. Hence, under near equilibrium condition i.e. ST being small, we observed that the cooling 
power is mainly dominated by interband processes by optical and SPP modes. Increasing the electronic temperature 
alleviates Pauli blocking, and allows for intraband processes to take place. As ST increases further, we observe that 
intraband cooling begins to dominate over the interband counterpart. The efficiency of energy exchange can be 
explained by the electron-phonon occupation number, quantified by the composite distribution function J-(Te,Tl) 
defined in Eq.|3] For inelastic processes, one can show that .F(Te,Tl) is independent of Tl when ST — > oo. This 
is as reflected in Fig.§for V OP ' K and V SPP ' H . On the other hand, for quasi-elastic acoustic phonons, the cooling 
power is proportional ST /Te instead. 

The results on doped graphene are shown in Fig.[2jb,d,f). Contrary to the intrinsic case, Pauli blocking promotes 
intraband electronic transitions over interband processes in doped graphene. In addition, "P" c ^ 'P" v , with larger 
cooling power for the majority carriers. At moderate doping of [i = 0.2 eV, their cooling power differs by more than 
an order of magnitude. The reduced electron-hole symmetry upon doping also leads to smaller interband cooling 
power. Quasi-elasticity of acoustic phonon scattering results in a phase space restriction in the scattering, with 
a Bloch-Gritneisen temperature determined by the doping 45, I60j. i.e. Tbg = 2hvskp /kp, in contrast to normal 
metals. This increase in phase space in conjunction with Pauli blocking greatly enhances the cooling power due to 
AP over the optical phonon baths. In fact, for moderate Te < 100 K, V AP dominates over all other mechanisms for 
cold graphene. 

The lattice temperature, Tl, also plays an important role in the competing cooling pathways. Fig. [3] compares 
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the fractional cooling powers V a jV T for intrinsic graphene, where V T — ^P 11 . To obtain a quantitative 
estimate, we include in-plane screening of the SPP scattering potential in graphene. The screening is incorporated 
through a standard procedure [rJPj |Mk,k'| — HMk.k'l / € 2d(<1, w). For simplicity, we employed the static screening 
dielectric function £2_o(q, 0), which in the long-wavelength limit assumes a simple form[62 e2£>(q, 0)~1 + q s / 'q, where 
q s —e 2 /2eoK J ^2?(e)de and T> is graphene density-of-states. / is the Fermi distribution function and is a function of 
the electronic temperature. 

We analyze the results in two non-equilibrium temperature limits, namely "near equilibrium" (Te — Tl = 10 K) and 
"far from equilibrium" (Te — Tl = 100 K) conditions. Fig. [3^ considers the condition of "near equilibrium" . At low 
Tl, AP dominates cooling. Increasing Tl populates the low-energy SPP mode, which begins to overtake the cooling 
power at a temperature of ~ 20 K. This transition temperature increases with doping e.g. is ~ 50 K at a doping of 
0.1 eV. The low energy SPP mode is overtaken by its high-energy mode at ~ 170 K. A downturn in the high-energy 
SPP cooling power is observed, due to larger screening at higher temperatures. Eventually, the optical phonons 
overtake the SPP for temperatures larger than 1000 K. Fig.[3j3 considers the condition of "far from equilibrium" . In 
this case, the SPP dominates the cooling power for all Tl, except at temperature >1000K where optical phonons 
begin to overtake it. 



B. Cooling dynamics 

We are interested in the role played by these various phonon baths on the cooling dynamics of photoexcited carriers, 
more specifically, the temporal evolution of Tg. The acoustic and optical phonon baths can each establish a different 
temperature upon interactions with the electrons, but processes such as anharmonic phonon-phonon scattering serve 
to thermalizc them on a picosecond time scale [51 IB5M65] . In this work, we shall assume a common lattice temperature 
Tl, but acknowledge that in experiments with ultrafast pump-probe, this will not hold true. On the other hand, 
under continuous light excitation, coupling to the heat sink via the supporting substrate substantially cools the 
lattice temperature to within a few degrees Kelvin of the ambient temperature T under usual photoexcitation 
conditions [66 . Typically, Te — Tl^> IK under low/moderate excitation power levels used in our studies. In this 
regard, the relatively small differences among the various phonon baths can be safely ignored. 

Hot carriers dynamics can be probed through optical measurements [TH51 167j. Following a pulsed light excitation, 
the temporal evolution of carrier relaxation, quantified by its electronic temperature Tg, can be measured using 
differential transmission spectroscopy. The electron dynamics are usually described by ATg cx exp(— t/re), and can 
be estimated with|68 , 



where C = d£ /dT^ is the electron specific heat and £ is the energy density of graphene. In this work, C is computed 
numerically. However, we note that for Tg <C fJi/ks, C increases linearly with T E , i.e. C « ^-V{p,)k 2 B T E - Having 
computed the total cooling power V T in Sec. Ill te can be obtained directly from Eq.|4j 



Fig. [4] plots the cooling time, te, for neutral graphene at Tl = 10 K. It is calculated for common substrates 
such as SiC>2, BN and non-polar substrate such as diamond. At very hot electron temperatures, i.e. Tb>500K, 
te is given by a relatively constant sub-picosecond cooling time. This is in agreement with experiments |22j . The 
constancy of te suggests the exp(— t/Ts) decay characteristics typical during the initial fast cooling process. As 
Te cools down, the cooling bottleneck due to AP begins to set in, leading to much slower cooling times. The 
transition temperature into this slow cooling regime varies with the choice of substrate as indicated in Fig. [4] This 
transition temperature is dictated by the lowest frequency SPP mode of the substrate. Unscreened results, which 
overestimate the SPP cooling power, yield much shorter cooling lifetimes than experimentally reported 22J. We also 
note that inclusion of disorder assisted cooling[68 , might enhance the decay rate, especially in the slow cooling regime. 

The optoelectronic response in graphene, photovoltaic [6"5H7T] or thermoelectric [2"2l [721 ITS] . is also a measure of the 
energy transport of these hot carriers. These experiments are usually performed under a continuous light illumination 
of an electrostatic junction. Their relative contribution depends on the electrostatic junction characteristics, doping, 
and even extrinsic factors such as electron-hole puddles[74]. Nevertheless, at steady state, the photovoltaic current 
is proportional to the photo-generated excess carrier density, Sn, via e<5n/Lt n £ where /i„ is the carrier mobility and 
£ the local electric field. The thermoelectric response, on the other hand, is proportional to the local elevated 
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temperature, ST=Te — Tl, via a(Si — S2)ST where Sx t 2 is the Seebeck coefficient of the two junction and a is the 
device conductivity. Here, we discuss estimates of Sn and ST. 

Under steady state condition, 

r° = J2v a +v M (5) 

a 

where V is the laser power absorbed by graphene, and V M is the heat dissipation via the metallic contacts, 
if any. In the absence of contacts, all heat dissipation is via the supporting substrate. At steady state, 
-pAP _|_ -pOP _|_ -pSPP (rp L _ t )k . Eq.JHjis then solved self-consistently in conjunction with charge conservation 
i.e. Sn — n e (TE,fi) — n e (7o,//o) =n h{TE, A 1 ) — nh{To, fJ*o)-> arriving at steady state values for Tg and /x. The 
photoexcited carrier density, Sn, and the elevated temperature, ST are plotted in Fig.[5|a-b) respectively, assuming 
typical experimental values of V = 1 x 10 7 W/m 2 and hq — 10MW/Km 2 . Both Sn and ST decrease with increasing 
ambient temperature T Q , due to more efficient cooling as phonon occupation increases. For SiC>2, Sn s» 10 10 cm -2 
and ST w 10 K under room temperature condition, of the same order typically seen in measurement [55] . Increasing 
doping increases the scattering phase space for intraband processes and lead to more efficient electronic cooling into 
the phonon baths, as shown in Fig. [2] As a result, both Sn and ST decrease with increasing doping. 

We also estimate the heat dissipation via contacts phenomenologically with V M Ri J(f — fo)(E — u)T> /tmcIE 
where fo is the distribution function before light excitation. First, we consider the simple ballistic limit where 
tm is just the device lifetime given by L/vf, where L is the length of the device and vf ~ 10 6 m/s is the Fermi 
velocity. Here, we assume a typical tm — 1 ps. We found that including V M only leads to few-fold decrease 
in the quantitative results presented in Fig. [5] (not shown). In the realistic case where the carrier transport is in 
the diffusive dominated regime, V M would be even smaller, by a factor of X/L, where A is the carrier's mean free path. 

Fig. [5j also suggests an order of magnitude enhancement in the optoelectronic response of graphene, by suppressing 
the SPP heat dissipation through a non-polar substrate, such as diamond-like carbon [75]. or by suspending graphene. 
In fact, the amount of heat transfer to the substrate via the electron coupling with the SPP can be quantified by an out- 
of-plane thermal conductance Kspp, defined as kspp = T' SPP /ST. This quantity sets the lower limit on the interfacial 
thermal conductance and it is plotted in Fig. § At room temperature, kspp ~ 1MW/Km 2 for undoped graphene, 
and can increases with doping to order of 10MW/Km 2 , see also Ref. [53]. For typical photocurrent experiments, L is 
typically 3> A, and transport is in the diffusive regime. Here, the in-plane electronic thermal conductivity, K e , can be 
estimated from the Wiedemann Franz relation. We found that our estimated value of kspp is significantly larger than 
K e /L 2 for typical experimental situations. This suggests that out-of-plane heat dissipation via SPP dominates over 
the in-plane electronic heat conduction. The former leads to an increased temperature of the graphene lattice. This 
result reconciles with recent experiment [66] , which reveals significant lattice heating upon laser excitation. From the 
experiment 66], we can estimate an out-of-plane thermal conductance of n exp = P°/ST ss 10MW/Km 2 . This value 
is consistent with our estimated kspp- In fact, klat alone is orders of magnitude smaller than the experiment as 
shown in Fig. [6] 

IV. CONCLUSIONS 

Our results point to the limiting role played by remote substrate phonons in the energy relaxation of hot 
photocarriers. In particular, we have shown that the steady state photoresponse in graphene is controlled by the 
inelastic scattering. The photovoltaic current is proportional to the photo-generated excess carrier density. The 
thermoelectric contribution, on the other hand, is proportional to the elevated electron temperature. Our results 
show that irrespective of the mechanism, the SPP phonons limits the overall strength of the photocurrent response on 
polar substrates. We predict that a choice of a non-polar substrate will lead to an order of magnitude enhancement 
in graphene photoresponse. Therefore, substrate engineering presents a promising route to efficient optoelectronic 
devices driven by hot carrier dynamics. 
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Appendix A: Acoustic phonon 

We consider first the energy exchange with the acoustic phonon (AP) bath. The total matrix element for electron- 
acoustic phonon scattering due to the two acoustic phonon modes, i.e. and Tta, is given bv|33[ [45] . 



\ M ^ = ^M (A1) 



where D ac is the acoustic deformation potential, taken to be 7.1 eV in our calulations, which is very similar to the 
recent ab-initio calculations of 6.8eV[76]. We note that the electron-phonon matrix element for these two acoustic 
modes have different angular dependencies with transition matrix elements |76l [77], which became negated after 
summing them [77]. vs is graphene effective sound velocity defined as [33] 2vg 2 = v^ A + v^ 2 Al where v$ = 17km/s, 
vla = 24km/s and vta = 14km/s. p m is graphene mass density taken to be 7.6 x 10~ 7 kg/m 2 . The acoustic phonon 
is then described by an effective Debye linear dispersion u q = v$q- Since vg <C Vp, hio q is typically much smaller than 
other energy scale in the problem. The acoustic phonon scattering is thus approximated to be elastic [IT)] i.e. k' k. 

The cooling power can then be written as, 

Vt P = #ffz r kdk [°° k'dk' f\o^^-8{k' -k-ujqvp 1 )^' -k^ik.k') 



(2n) 2 hJ J J 2p 

7-1 2 fc POO p27T 

9s9v D z ac h ■ 



rOO r-ZTT 

/ / k 2 q 2 —j r (k,k)dkd9 (A2) 
Jo Jo v F 



(2n) 2 h 2p m v s Jo Jo 
Under the assumption Hui q <C Te,Tl, we have 

JXM)«(i-A)A^^ (A3) 

Making use of the relations q 2 ~ 2fc 2 (l — cos9), we then obtain a simplified form for the cooling power, 

Vt P * 9 * 9 « D «' Te - Tl (°° k\l - f k )f k dk (A4) 



2irp m vs T E jo 

Contributions from interband processes, V^ p vc , are forbidden due to energy-momentum conservations 



Appendix B: Optical phonons 

Next, we discuss energy exchanges with high energy dispersionless optical phonon (OP) modes i.e. TlO: ^to an d 
K-pq. We consider first the electron-phonon coupling of long- wavelength optical phonon modes, T^o and Tto- Their 
sum is expressed as [161 ITT] . 



A I, 



OPT 
k,k' 



nn 2 

lllJ OPT 

2p m io 



(Bl) 



where Dop.t — iV2g/2 m lleVA 1 is the optical-phonon deformation potential with a coupling constant of 
g = 5.3eVA _1 , and Huj = 197 meV. We note that the electron-phonon matrix element for these two optical modes 
have different angular dependencies with transition matrix elementsjlSl 27] i.e. 1 ± ss' cos(#k — 8^) where s = ±1 
denotes conduction/ valence bands, which again became negated after summing them. 



We consider first the intraband cooling power, written as, 

OPT _ 9s9v 



V. 



r°° r 2n D 2 h 
. . kdk I k'dk' / d9 OP ' T S(k' -k-LJovZ^k' -k)T(k,k') 
(2iryhJo Jo Jo 2pm.ua 

9s9vDqpy 
9s9vD 



k{k + ojov F 1 )j r (k, k + uiov F 1 )dk 



2 

OPT 



4irp m v F 



exp( 



hujo 

kpTh 



-exp(— — ) 
KbTe 



k(k + wou/Xl - fk)f k+LJnV -idk 



(B2) 
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In similar fashion, the interband cooling power is written as 

g s gvD 2 



V, 



OPT 



'OPT 



4irp m v F 

where H v is the Heaviside function. 



kb-Ll kb±e 



fc(w(V - k)H v [oo Q v' F l - k](l - f k ),f k _-!dk (B3) 



For zone edge phonon modes, only the transverse K.pQ contributes to carrier scattering, and the matrix element 



is 



2 _ hD 2 OP K i - ss 1 cos e 

2p m UJ o 2 



(B4) 



where Dop,k = 3<? ~ 16eVA 1 and Hujq — 157 meV. The cooling power is similar to the T phonons case, except a 
factor of \ smaller due to the angular dependence. 



Appendix C: Surface polar phonons 



The surface polar phonons coupling is given by 26, 32, 33 



SPP\ 



-F 



2 exp(— 2qzo) 1 + ss'cos9 



(CI) 



where eo is the free space permittivity and zq is the separation between graphene and the substrate. The magnitude 
of the polarization field is given by the Frohlich coupling parameter, F 2 . In common SiC>2 dielectrics, there are 
two dominant surface optical phonon modes having energies fvjj\ — 58.9 meV and hu)2 = 156.4 meV, with Frohlich 
coupling Fl = 0.237 meV and F 2 2 = 1.612 meV respectively [53] . 

We consider first the intraband cooling power, written as, 



,SPP 



9s9v 

(2nynj 



kdk 



k'dk' 



(10 



ne 2 2 exp(— 2qz ) 1 + ss'cos9 , 



S(k' -k- UjVp 1 )^' - k)T(k, k') 



g s g v uij we 2 
{27r) 2 hv F e 3 



k(k + ujjV F 1 )J r (k, k + ujjVp 1 ) 



exp(— 2qz ) 1 + cos9 



dOdk 



(C2) 



The phonon momentum q has the constraint q 2 = k' 2 + k 2 — 2k'kcos9. Under typical conditions, the factor 
exp(— 2qz(j) rj 1. Linearizing exp(— 2qzo), the intraband cooling power then becomes, 



V SP P = !!-!!■ ~~' 2 '''>y_ 
(2ir) 2 hvpeQ 

where k' = k + ujQVp 1 and 



kb-Ll k B J-E 



kk'{l-f k )f k ,G(k,k')exp 



-zq 



Q(k,k') 



- dk (C3) 



Q(k,k') 



1 + COS0 



o 2^/k' 2 + k 2 ~ 2fc'fccos<9 



.d9 



k + k' 
kk' 



Ik 



2Vkk' 
k + k' 



Ie 



2Vkk' 
k + k' 



(C4) 



where Ik,e are the complete elliptic integrals of first and second kind. In similar fashion, the interband cooling power 
is written as, 



, S pp _ <i>nr~j-> 2 r; y 



kB-LL ksJ-E 



{2iz) 2 hvF£o 

jf kk'H v [k'](l - h)f- k ,Q{k, -fc')exp f e( ~^ fc0 



dk 



(C5) 



where kl =ujqv f x 
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FIG. 1: (Color online) Cartoon illustrating the typical cooling pathways of hot electrons produced by continuous photoexci- 
tation of graphene with detail descriptions in the main text. Heat can also be dissipated through metallic contacts attached to 
graphene (not shown), as discussed in Sec. Ill B 
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FIG. 2: (Color online) (a-b) Electron cooling power due to acoustic phonons bath, V p , as function of electron temperature 
Te for different lattice temperature Tl calculated for neutral and doped graphene respectively. Intraband (cc, vv) and interband 
(cv, vc) processes are indicated, (c-d) Similar, except for K optical phonons modes, r p OP ' K and (e-f) for high energy unscreened 
SPP mode, -p SPP ' H . r optical and the low energy SPP phonons show similar characteristic (not shown). 
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FIG. 3: (Color online) (a) Fractional cooling power V a /V T where V T = E Q calculated at 1 'near equilibrium"' condition 
of Te — Tl = 10 K for neutral graphene, including 2D screening £20(9) described in text. For SPP and OP, the dashed 
(solid) line represents the low (high) energy mode, (b) Same as (a), except calculated for "far from equilibrium"' condition of 
Te-Tl = 100K. 
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FIG. 4: (Color online) Carrier's cooling time te after photoexcitation plotted as function of Te calculated for cold neutral 
graphene. Various substrates are considered, namely Si02, BN and non-polar, calculated for screened and unscreened SPP 
scattering potentials. Dotted line distinguishing fast/slow cooling is rather arbitrary, and serve only as guide to the eye. 
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FIG. 5: (a) (Color online) Steady-state excess carrier density, 5n, upon continuous photoexcitation as function of the ambient 
temperature for neutral graphene. Various substrates are considered, namely Si02, BN and non-polar, calculated for screened 
and unscreened SPP scattering potentials, (b) Elevated temperatures, Te — Tl, calculated for same conditions in (a). All 
calculations assumed P° = 1 x 10 7 W/m 2 and k = 10MW/Km 2 . 
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FIG. 6: (Color online) Out-of-plane thermal conductance, kspp, denned as kspp = V /ST calculated for Si02 for different 
conditions such as (i) screened and unscreened SPP scattering potentials and (ii) different doping fj, (all curves are for zero 
doping unless stated otherwise), klat for undoped graphene on a non-polar substrate is plotted as reference. 



